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The sub-stellar birth rate from UKIDSS 
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ABSTRACT 

We present a new sample of mid L to mid T dwarfs with effective temperatures of 
1100 to 1700 K selected from the UKIDSS Large Area Survey and confirmed with 
infrared spectra from X-Shooter/VLT. This effective temperature range is especially 
sensitive to the formation history of Galactic brown dwarfs and allows us to constrain 
the form of the sub-stellar birth rate, with sensitivity to differentiate between a flat 
(stellar like) birth rate, and an exponentially declining form. We present the discovery 
of 63 new L and T dwarfs from the UKIDSS LAS DR7, including the identification of 
12 likely unresolved binaries, which form the first complete sub-set from our program, 
covering 495 sq degrees of sky, complete to J=18.1. We compare our results for this 
sub-sample with simulations of differing birth rates for objects of mass 0.10-0.03M Q 
and ages l-10Gyrs. We find that the more extreme birth rates (e.g. a halo type form) 
can likely be excluded as the true form of the birth rate. In addition we find that 
although there is substantial scatter we find a preference for a mass function, with a 
power-law index, a in the range — 1 < a < that is consistent (within the errors) 
with the studies of late T dwarfs. 

Key words: low mass stars, brown dwarfs. 



1 INTRODUCTION 

The distribution of star formation with mass and time 
are key pieces of observational evidence for understanding 
star formation in the galaxy. The former is d escribed by 
the initial mass function (IMF; ISalpetenll995l ). which can 
be described as a power- law of the form ip(M) oc M~ a , 
with a — 2.35, and has been determined across the stellar 
mass regime by measuring the luminosity function for 
a population of stars, and applying a mass-luminosity 
relation, which should account for metallicity variations. 
Since brown dwarfs never reach the main sequence, this 
determination is complicated in the sub-stellar regime 
by the lack of a unique mass-luminosity relationship. 
Instead t he T e s and lumin osity are dependent on mass 
and age l|Allard et al.l I1997T ). This means that the lumi- 
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nosity function and T e g distributions of field brown dwarfs 
depend not only o n the mass fun ction, but also on their 
formation history i|Chabrierl 2002 V Indeed, extending the 
Salpeter mass function to sub-stellar objects one would 
expect many more brown dwarfs than stars, which is 
not seen by observations of late M and L dwarfs (e.g. 
iReid. Gizis fc Hawlevl 12051 ) . In addition it is not totally 
apparent that the full field population of brown dwarfs is 
populated by objects that formed through a cloud frag- 
mentation process. There are a number of ways that brown 
dwarfs could form that are different to the canonical ways 
in which stars are thought to form. There are at least four 
different formation mechanisms that have been suggested, 
such as formation through the "ejection" of pre-stellar 



20011; iBate. Bonnell fc Bromm 



cores ( Rcipurt h fc Clark 
| 2002t iDelgado-Donate. Clarke fc Bate! 120031 ; ISterzik et al 
120031 ) o r through "turbulen c e", o r turbulent fragmen- 



tation i|Padoan fc Nordlundl |2002| ; iPadoan fc Nordlundl 
12004 ). Other theories include that of "disc fragmenta- 
tion", forming sub-stellar cores from an initially massive 
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pre-st ellar core via fragm e ntation of a large circumstellar 
disk dBoffin et all Il998l; iBate, Bonnell fc Bromml 120031; 



Whitworth fc GoodwinTl2005l : I Whitworth fc Stamatellod 
20061 ) and "photo-erosion", where sub-stellar objects form 



in the presence of a higher mass star embedded in a HII 
region ijWhitworth fc Zinneckerl 120041 ') . As the formation 
mechanisms could indeed be different for sub-stellar objects 
it is thus important to define the mass function and for- 
mation history in the sub-stellar regime if we wish to fully 
understand their contribution to the Galactic population. 

Young clusters have been the target of many stud- 
ies seeking to measure the sub-stellar IMF since their 
known ages and metallicities allow the use of a mass- 
luminosity rela t ion based on a single coeval age (e.g . 
Caballero et~aTl 120091; lOliveira. Jeffries fc van Loonl 120091 ; 



Luhman et all 120071 ; lLodieu et"ai1 1200781 ). Although these 



clusters allow a relatively direct measurement of the sub- 
stellar IMF, they also introduce their own problems, since 
the initial conditions and accretion histories of individ- 
ual objects introduce uncertainties regarding th e ages, 
and hence masses, of such young objects (e.g. iBaraffd 
l2010l ). As it is difficult to determine the age of field 
brown dwarfs, unless they have fiducial constraints on their 
age a s binaries (e.g. Bur ningham ct al. 2009; Zhang ~et all 
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20101; iBurningham et al.l l2010al ; iDav- Jones et all 1201 ll; 



Murr ay et al. _ 2011 ) or as members of moving groups 



l2010l ; IClarke et al.ll2010l ). estimating the 



mass function of field brown dwarfs requires a knowledge 
of their formation history. This is often ass umed to be the 
same as that for stars (constant with time; iMiller fc Scald 
1979), but is unconstrained in the sub-stellar regime. The 
first atte mpt to mea s ure th e sub-stellar mass function was 
made by iReid et "all (|l999T l from preliminary results from 
2MASS data based on only 17 1 dwarfs. More recently, 
with the discovery of late T dwarfs several other groups 
have made measurements of the sub-stellar mass function in 
the disc. These have all generall y been with smal l sample 
sizes or cover only 1 dwarfs (e.g. Cruz et al. 2007 ) or only 
T dwarfs (e.g. [Burningham et all 2010bl ; IKirkpatrick et all 



120121 ; iMetchev et al.ll2008l ). Those that have considered the 
full temp erature regime acr oss the L and T dwarf spectral 
types (e.g lRevle etaTll2010l ) suffer from large associated er- 
rors and large bin sizes in order to get large enough sampling. 

In order to characterise the form of the sub-stellar 
formation history a large sample of brown dwarfs is re- 
quired. With modern large -scale near and mid in frared sur- 
yeys, such as th e 2MA SS dSkrutskie et al.ll2006l). UKIDSS 
ilLawrence et all |2007|). VISTA (lEmerson fc Sutherland! 
120031 ) and WISE i|Wright et aT1l2010l), which have identi- 
fied large numbers o f brown dwarfs (IKirkpatrick et al. 2000l; 
Hawlev et all |2002|; IPmfield et alJl200St (Burningham et all 
2010b ; Kirkpat rick et alll201l] . ICushing et al. 201l]) it is now 



possible to provide the necessary sample of such objects. 
This paper outlines our efforts to use the UKIDSS large 
Area Survey to empirically constrain the Galactic brown 
dwarf formation history. We discuss past and present sim- 
ulations of the formation history in 321 In [J3]the selection 
of our sample of 1 and T dwarfs from the UKIDSS IAS, g] 
gives details of the observations and data reduction. In fj5] 
we present the spectroscopy and spectral types of our sam- 
ple and investigate potential unresolved binarity. In ij6] we 
compare our first observations with simulations and look at 



constraints that can be placed on the formation history of 
Galactic brown dwarfs. Finally we summarize our findings 
in § El 



2 SIMULATIONS OF THE FORMATION 
HISTORY 

Simulations of the effect of how the birth rate affects the lu- 
minosity function have been performed by several authors. 
Based on model data an d direct compar isons with DENIS 
and 2MASS observations. IChabrierl (|2002j) made simulations 
using two different initial mass functions and birth rates. 
Namely a flat, or constant birth rat e, which is the s i mples t 
form and supported by the work of lMiller fc Scald (|l979l ). 
who suggest that the birth rate does not depend strongly 
on the gas density, and is approximately consistent across 
the Galactic disc. They also consider an exponential form 
where the formation decreases with time. These are con- 
sidered with initial mass functions of a power-law form de- 
rived bv lChabrieil(|200ll ). and a log normal and exponential 
form which essentially give the sam e result wh e n con sider- 
ing the effects from the birth rate. iBurgasseil l|2004h con- 
sidered a wider range of birth rates in his monte-carlo sim- 
ulations considering, in addition to a flat and exponential 
forms, an "empirical" birth rate, which i s the s ame as that 
measured for stars bv lRocha-Pinto et all (|2000| ). which rep- 
resents 'bursts' of formation at peak internavls of 0-1, 2-5 
and 7-9Gyrs. This formation history scenario is also sup- 
ported by the more re c ent w ork o f the stellar f ormation his- 
tory bvlCignoni et al l |2006h and lWvsd (|2008l ). In addition 
iBurgassei 1 20041 ) also considers a "cluster" birth rate which 
assumes a flat, but stochastic (i.e in a number of clusters) 
formation, which produces a similar result to a flat forma- 
tion scenario. Finally he considers a "halo" type birth rate, 
that includes formation within a 1 Gyr burst, 9 Gyr in the 
past, in an attempt to explain a num ber of subdwarf brown 
dwarfs that have b een identified (e.g. iBurgasser et al .1120031 : 
llodieu et alj|2O10l ). This scenario gives a radically different 
T e ff distribution for L and T dwarfs compared to the other 
scenarios, and seems unlikely since we are now seeing a larger 
number o f L dwarfs identifie d in very young clusters (e.g. 

Quanz et ai1l2010l . Cameleo n; lluhm an 200^ , Ser- 

;IChauvin et al.l l2 



Taur us; ?,^, _ _ . . _ . 

pens; llodieu et al.ll2002l TWA Hydra; Ichauvin et alj|2004 



Upper Sco; lLodieu et aTll200Sl . iLafreniere et al.ll2008l ) . 

More recen t simulations performed by 

iDeacon fc Hamblvl (|2006l ) looked more specifically at 
1 and T dwarfs from the UKIDSS LAS. They produced 
simulations that take into account several IMFs including a 
flat, log normal and different power laws (a =+1.0, 0, and 
-1.0), combined with different exponential forms of the birth 
rate, similarly to those described above. These simulations 
also included the effect of Galactic disc heating, which 
had not been included in previous simulations of the birth 
rate. A histogram of these simulations (a log normal form 
of the IMF with different birth rates) i s shown in Fig. [fl 
These simulations are similar to those of lAllen et all (| 19951 ) 
(see their fig. 2). with the main differences arising from 
the differen c es in the normalisation of the space density. 
I Allen et all (1 19951) use s 0.35 stars pc 3 Mq according 
izid lli 



to 



Reid fc Gizid ll 19971). where as the simulations based 



on those of IDeacon fc Hamblvl (|2006l ) use 0.0024 stars 
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pc -3 Mo according to lDeacon. Nelemans fc Hamblvl {2008). 
In addition changes also arise from the different values used 
for their birth rates, and as such are similar and show the 

same trends but are not directly comparable. 

It can be clearly seen in the iDeacon fc Hamblvl l|2006h 
simulations that the sub 1000K region is extremely sensitive 
to the IMF, but relatively insensitive to the birth rate, while 
the 1100-1500K, corresponding to the mid-L to mid-T spec- 
tral range is most sensitive to differing birth rates. While 
several hundred brown dwarfs have been identified in large 
area surveys such as 2MASS and SDSS, they were more sen- 
sitive to the detection of L dwarfs and produced only a few 
tens of early T dwarfs. As such they could not provide the 
population needed to study and constrain the birth rate. 
The UKIDSS LAS probes to greater depth across the L and 
T dwarf spectral types and can provide a statistically robust 
sample spanning the mid-L to mid-T region, which is most 
sensitive to the effects of the form of the birth rate. We thus 
select a sample of mid L - mid T dwarfs from the UKIDSS 
LAS in order to compare the s pace density with that of late 
T's l|Burningham et al]|2010ah and late Ms and earlier Ls to 
measure the birth rate. 

3 SAMPLE SELECTION 

We select our sample from DR7 of the UKIDSS LAS survey, 
selecting objects with declination ^ 20 degrees, J < 18.1 an d 
Y - J > 0.8 following the criteria of iHewett et alJ l|2006h . 
who demonstrate that using a Y — J > 0.8 colour criterion is 
largely free of M dwarfs while including all L dwarfs £16.11 We 
also used the following quality flags to ensure that selected 
objects are point sources and are not likely cross talk effects, 
or sit at the edge of the detector, such that the following 
quality flags were used: 

(priOrSec= OR priOrSec = frameSetID) 
yppErrBits < 256 
jippErrBits < 256 
hppErrBits < 256 
—3.0 < mergedClass < —0.5 
—3.0 < mergedClassStat < 3.0 
yEll < 0.45 
jilEll < 0.45 

This list of candidates was then cross- matched against SDSS 
DR7 to identify objects with optical counterparts. Using a 
matching radius of 4 arcsec we selected objects with opti- 
cal cou nterparts according to the criteria below, based on 
those of I Schmidt et al.l (|2010l ). who provide colours from an 
unbiased spectroscopically complete sample. Since our NIR 
colour selection effectively removes all contaminant field M 
dwarfs this allowed us to be more liberal with our redder 
sources (J — K > 1.0) in terms of our z — J colour selec- 
tion, and allows for larger uncertainties in the i' band as 
we are probing the faint end of SDSS. We also required Y 
and J photometry to have errors not greater than 3a and 
detections in H ^ 14.5, such that this search space would 
not have been probed by 2MASS. In addition we also con- 
sider K band non-detections, if their z — J colour passes our 
following criteria: 

J < 18.1 

Y - J ;s 0.8 



z - J ^ 2.4 and ( J - K > 1.0 or no K detection) OR 
z — J ^ 2.9 and (J — K < 1.0 or no K detection) 
We then removed objects with the following: 
i - z < 2.0 and a(i - z) < 0.35 
i - J < 4.7 and a(i - j) < 0.2 
z - K < 3.5 and J - K > 1.0 and a{z - K ) < 0.2 

As mid-T dwarfs typically have z' — J > 3.0 
fe.g lPinfield et all [2008) some objects will be too faint 
for detection in SDSS, we thus include these SDSS non- 
detections. All objects were then visually inspected to re- 
move any possible miss matches or cross-talk. We also cross 
matched our sample with known L and T dwarfs in dwar- 
farchives and retain those that have spectral types of >L4 
(to be consistent with our completeness, see i)6.1[) . giving 
a total sample size of 324 from 2000deg 2 of UKIDSS LAS 
(DR7) sky. In this work we consider a sub-sample from a 
smaller are of 495 sq degrees of sky in the RA and DEC 
range, RA=22 to 4 hr and DEC=-2 to 16 deg. This includes 
76 L and T dwarfs, 13 of which were previously identified 
(and spectrally typed as L4-T4), 5 of these have been re- 
observed and 63 new L and T dwarfs. We show photometric 
information of this sub-sample in Fig. [2] and Table. [1] which 
are discussed and analysed further in the following sections. 



4 OBSERVATIONS & DATA REDUCTION 

Spectroscopic observations of candidates from our sub- 
sample, were obtained with X-shooter on the Very Large 
Telescope during Nov 27-31 2010, Feb 22-26, June 5-9 
2011 and Sep 18-22 2011, under the ESO programs 086.C- 
0450(A/B) and 087.C-0639(A/B). We used the echelle slit 
mode, which covers the wavelength range 300-2500nm. This 
is split into three separate arms, the UVB (300-550nm), VIS 
(550-1000nm) and NIR (1000-2500nm). Using slit widths of 
1.0 arcsec for the UVB arm and 0.9 arcsec for the VIS and 
NIR arms, we took four individual integrations in an ABBA 
pattern (See Table. [2] for integration times). We note that 
we do not detect any significant flux in the UVB wavelength 
range for members of our sample and as such we do not show 
the details for this arm. We took telluric standard stars af- 
ter every second target, which were paired together so they 
were in roughly the same airmass. Sky flats and arc frames 
were also taken at the beginning of every night. 

The data were reduced using the ESO X-shooter 
pipeline (version 1.3.7). The pipeline removes non-linear pix- 
els, subtracts the bias (in the VIS arm) or dark frames (in 
the NIR arm) and divides the raw frames by flat fields. Im- 
ages are pair-wise subtracted to remove sky background. 
The pipeline then extracts and merges the different orders 
in each arm, rectifying them using a multi-pinhole arc lamp 
(taken during the day-time calibration) and correcting for 
the flexure of the instrument using single-pinhole arc lamps 
(taken at night, one for each object observed). Telluric stars 
are reduced in the same way, except that sky subtraction 
is done by fitting the background (as tellurics are not ob- 
served in nodding mode). The spectra were telluric corrected 
and flux calibrated using IDL routines, following a standard 
procedure: first the telluric spectrum is cleared of HI ab- 
sorption lines (by interpolating over them) and scaled to 
match the measured magnitudes; then is divided by a black- 
body curve for the appropriate temperature, to obtain the 
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Table 1. Photometry of sub-sample. 
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:06: 


04 


T"Jd.^±y.oy 


1 7 


Q9 

yz 


1 


.06 


^0 


1 11 


2 


,18 


1 


.92 


SDSS T09074-0000 a 

OUOO dU^UI tUUUU 




02 


:07, 


:42 


-UOO-OO-^fi 

TvW ■ UU. Ok) 


16 


73 


1 


.29 


-0.07 


0.01 




_ 






TTT, AS T0947 0107 


RR T ,T4S 


02 


:47: 


03 


-01:07:01 


17. 


.77 


1 


.43 


0.94 


1.77 




_ 






TTT, AS T09 c i c i4-0fi1 H h 


RRT,T^0 


02 


:55: 


:45 


4-Ofi-l f\-^f\ 
T"U ■ -LU. UU 


17. 


99 


1 


.16 


-0.68 






_ 






TTT AS T09 c ;Q-U0 c ;4Q 


RRT T^l 


02 


:59: 


Al 


4-0 ^- 40-^^ 


1 8 

lo, 


09 
uz 


1 


.25 


ft 1 } 

U.OJ 


1 ^4 
1.04 




_ 






TTT, AS T0*31 44-04^ 


RRT,T^9 


03 


:14; 


52 


4-04- c i^-4fi 

U^i . U O . ^4: u 


17. 


30 


1 


.29 


0.91 


1.71 


3 


.06 


2 


.00 


TTT, AS T0^904-0(S1 7 


RRT,T^fi 

±_>±\,J_J _L OKJ 


03 


:20 


00 


4-0(S-1 7-41 


17, 


79 


1 


.46 


0.75 


1.39 




_ 






TTT, AS T0^904-0fi1 8 


RRT,T^7 

J_> 1 V J_J _L O 1 


03 


:20: 


■42 


4-0fi-1 R-^7 


18 


06 


1 


.21 


0.58 


1.15 


2 


.43 


1 


.93 


TTT AS T0^91 4-0^4^ 

Ui-JiT-D J UOZ !_-(-- U 040 


RRT T^R 

D IV Lj _L OO 


03 


:21: 


:43 


4-0^-4^-94 


1 7 


OO 


1 


.22 


7 1 } 


1 ^7 


2 


.74 


2 


.34 


TTT, AS T0^9^4-0(S1 ^ 


RR T ,TfiO 

J_> JTLJ-i -L \jyj 


03 


:23 


:54 


4-0fi-1 ^-^9 


17 


64 


1 


.37 


0.65 


1.33 




_ 






SDSS T0^9 c i4-049 c i/ 




03 


:25: 


53 


4-04-9^-40 

TU4.ZiO.4U 


16 


02 


1 


.10 


-0.22 


-0.43 




_ 






TTT, AS T0^'^04-0 r i^fi 


RRT,Tfi9 

J_> JTLJ-i J- 


03 


:30: 


06 


TuJ •Okj.OO 


17 


95 


1 


.56 


1.10 


2.00 




_ 






TTT AS T0^^04-049fi 


RRT Tfi4 

O IV Lj J- U4: 


03 


:30: 


:37 


4-04-9(^-^8 


1 7 


90 
zy 


1 


.33 


£^ 

U.OJ 


1 ^4 
1.04 


2 


.69 


2 


.16 


TTT-AS T0^41 4-042^ 


J_J J. VU 1 uu 


03 


■Al: 


■50 


4_04-2'>-2 c i 


16. 


.85 


1 


.44 


0.90 


1.65 


2 


.93 


2 


.12 


TTT, AS T91 ^74-00^8 


RRT^^O^ 

JJIVU X ijKJU 


21 


:57, 


00 


4-00- ^fi-1 ^ 
T«U . UU. -LU 


17 


85 


1 


.41 


0.98 


1.75 


2 


.77 


1 


.66 


TTT, AS T91 ^Q-I-OO^ 1 } 


RRT,T^0fi 

JJJ.V.J.J J- ouu 


21 


:59: 


20 


4-00-^^-1 

T«U . OO. -LU 


17. 


73 


1 


.36 


0.74 


1.37 


2 


.76 


2 


.51 


TTT AS T990Q 00^^ 


RRT T^07 

D IV 1j J. OkJ I 


22 


:09: 


:17 


00- ^3-00 

-uu.oo.uu 


1 8 

lo. 


1 


1 


.34 


7 A 

U. / 4 


1 ^7 


2 


,18 


2 


.31 


ULAS J2229+0102 


BRLT311 


22 


:29; 


:58 


+01:02:17 


17. 


.88 


1 


.22 


0.39 


0.67 










ULAS J2233+0022 


BRLT312 


22 


:33: 


:48 


+00:22:14 


18. 


07 


1 


.05 


0.71 


1.43 


3 


.45 





.42 


ULAS J2236+0111 


BRLT313 


22 


:36: 


37 


+01:11:32 


17. 


11 


1 


.34 


0.87 


1.64 


3 


.04 


2 


.03 


ULAS J2237+0716 


BRLT314 


22 


:37: 


:57 


+07:16:57 


17. 


49 


1 


.38 


1.04 


1.84 


2 


.97 


2 


.73 


ULAS J2240+0008 


BRLT315 


22 


:40: 


52 


+00:08:22 


17. 


.82 


1 


.46 


0.70 


1.24 


2 


.44 


2 


.03 


ULAS J2249+0715 


BRLT316 


22 


:49 


23 


+07:15:28 


18. 


09 


1 


.55 


0.55 


1.23 










ULAS J2250+0808 


BRLT317 


22 


:50: 


16 


+08:08:22 


15 


50 


1 


.17 


0.46 


0.99 


2 


.74 


2 


.11 


ULAS J2251-0007 


BRLT318 


22 


:51: 


:15 


-00:07:24 


17. 


95 


1 


.26 


0.60 


1.46 


3 


.20 





.97 


ULAS J2256+0724 


BRLT320 


22 


:56: 


31 


+07:24:39 


17. 


94 


1 


.48 


0.68 


1.21 


2 


.71 


2 


.45 


ULAS J2302+0700 


BRLT321 


23 


:02: 


03 


+07:00:39 


17. 


62 


1 


.33 


0.25 


0.11 










ULAS J2303+0058 


BRLT322 


23 


:03: 


:59 


+00:58:07 


17. 


82 


1 


.21 


0.83 


1.67 


2 


.86 


2 


.59 


ULAS J2304+1301 


BRLT323 


23 


:04: 


25 


+ 13:01:11 


16. 


69 


1 


.31 


0.77 


1.49 


2 


.78 


2 


.06 


ULAS J2304+0804 


BRLT325 


23 


:04; 


34 


+08:04:01 


17. 


89 


1 


.23 


0.41 


0.67 
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Figure 1. Monte-Carlo simulations of differi ng IMFs, where i/)(M) oc M Q for values of a = +1.0,0.0 and —1.0 and birth rates 
/3 = 0.0,0.2,0.5, based on the simulations from lDeacon fc"H amblv ( 2 0061 ). 



Table HI continued. Photometry of sub-sample. 



Name 




ID 


RA 


DEC 


J 1 


Y-J 1 


J-H 1 


J-K 1 


z-J 1 - 2 


i-z 2 


ULAS 


J2306+1302 c 




23:06:01 


+13:02:25 


17.57 


1 


39 


-0.43 


-0.46 






ULAS 


J2312+0006 


BRLT328 


23:12:37 


+00:06:02 


17.65 


1 


30 


0.60 


1.25 


2.62 


2.07 


ULAS 


J2316+0100 


BRLT330 


23:16:46 


+01:00:13 


17.95 


1 


15 


0.69 


1.25 


3.13 


1.94 


ULAS 


J2320+1448 9 




23:20:35 


+14:48:29 


16.79 


1 


35 


-0.35 


-0.61 






ULAS 


J2321-0045 


BRLT331 


23:21:23 


-00:45:57 


18.00 


1 


10 


0.40 


0.88 


2.52 


2.44 


ULAS 


J2321+13549 




23:21:23 


+13:54:54 


16.72 


1 


20 


-0.43 


-0.44 






ULAS 


J2323+0005 


BRLT332 


23:23:00 


+00:05:42 


18.01 


1 


15 


0.74 


1.15 


2.53 


2.17 


ULAS 


J2323+0719 


BRLT333 


23:23:15 


+07:19:31 


17.30 


1 


20 


0.75 


1.10 


3.05 


3.38 


ULAS 


J2327+1517 


BRLT334 


23:27:16 


+15:17:30 


16.20 


1 


34 


0.85 


1.52 


2.97 


2.13 


ULAS 


J2327+0102 


BRLT335 


23:27:32 


+01:02:53 


18.07 


1 


19 


0.83 


1.46 






ULAS 


J2328+1345 c 




23:28:02 


+13:45:44 


17.75 


1 


26 


-0.42 


-0.54 






ULAS 


J2330+1403 


BRLT338 


23:30:02 


+14:03:30 


17.37 


1 


23 


0.57 


1.26 






ULAS 


J2347-0110 


BRLT343 


23:47:17 


-01:10:09 


17.57 


1 


25 


0.85 


1.67 


2.70 


2.35 


ULAS 


J2356+0754 


BRLT344 


23:56:18 


+07:54:20 


18.09 


1 


51 


1.10 


1.87 







I : Colours and photometry from UKIDSS (Hcwctt ct al. 2006), 2 : Colours and photometry from SDSS (Fukugita ct al. 1996). 
Original discovery paper: "Geballe et al. (2002) . 6 Hawlev et al. (2002) . c Burningham et al. (2010a), d Knapp et al. (2004) . 



e Lodieu et al. (2007a) . ^Chiu et al. (2006) . s_Scholz (2010) . h Burningham et al., in prep. 



instrument+atmosphere response curve; finally the target 
spectra is multiplied by the response curve obtained to flux 
calibrate it. The arms (VIS and NIR) were then merged by 
matching the flux level in the overlapping regions between 
them. The flux calibration was checked by determining the 
target's synthetic MKO YJHK magnitudes, that were com- 
pared to those obtained in the ULAS. Finally, each spectrum 
was visually inspected to check for possible problems dur- 
ing the extraction or merging stage. The spectra were then 
binned (in the A direction) by 40 times to produce an aver- 
age S/N= 30 for resolution R= 880 and 510 in the VIS and 
NIR arms, respectively. 



5 SPECTROSCOPY 
5.1 Spectral Types 

Spectral types were determined via spectral fitting with 
template objects. Template spectra were obtained from the 
SpeX Prism Spectral Library 0, which are shown in Ta- 
ble. [3] Our observed spectra were then re-sampled to the 
same resolution as the SpeX data (see Table. [3} and spec- 
tral types were determined from the best chi-squared fit 



1 http://pono.ucsd.edu/ adam/browndwarfs/spexprism/ 
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Table 2. Optical and near-infrared spectroscopic observations of sub-sample. 



Name 


UT date 


VIS t i7lt (s) 


VIS S/N* 


NIR t int (s) 


NIR S/N* 


BRLT1 




1600 


4.54 


1960 


29.7 


BRLT2 


90 ^lort 901 1 


2000 


4.72 


2360 


18.8 


BRLT3 

!_} 1 LJ-J X o 


28-Nnv-201 


1200 


5.21 


1560 


32.8 


RRT Tfi 

J_> 1 1_Xj _L U 


1 8 ^on 901 1 


2000 


5.32 


2360 


25.2 


BRLT7 


9(1 ^on 901 1 


1600 


7.61 


1960 


20.9 


RRT TR 


91 CJnr» 901 1 


1 fiOO 
1UUU 


^18 
o. xo 


1 QfiO 
xyuu 


31 

Ox 


RRT TQ 

j_> ± _i_ y 


1 Q ^on 901 1 

X y Ot;JU— X 


1600 


8.01 


1960 


30.7 


RRT T1 

J_> 1 X, J_j 1 1U 


9Q Mnv 901 


800 


2.39 


1160 


26.5 


RRT T1 9 

J_> ± V.-LJ X X Zr 


91 ^pn 901 1 


2000 


4.3 


2360 


18.7 


RRT T1 4 

J_> J. \,_L; _L J. 1 * 


30 Mnv 901 


1600 


7.56 


1960 


29.8 


RRT T1 

JJlilj X XO 


90 ^on 901 1 


1600 


2.4 


1960 


21.6 


RRT.Tl fi 

J_> ± i_Xj 1 1U 


98 Mnv 901 

ZO 1 N U V ZiUlU 


1400 


3.66 


1760 


22.8 


RRT T1 8 

UIUj X _LO 


30 Nnv 901 


1600 


5.79 


1960 


29 


RRT T90 


91 CJnr» 901 1 


9000 
ZUUU 


5 12 


93fi0 

zouu 


90 7 


RRT T91 


1 Q Son 901 1 

X 3 JCU illl X 


1200 


3.73 


800 


29.5 


RRT T99 

J_> ± l_Xj 1 Zi^ 


91 Son 901 1 

Z 1 UCU iUl X 


1600 


5.84 


1960 


17.7 


RRT T94 


91 Son 901 1 


1600 


4.17 


1960 


23.9 


RRT T9fi 

J_>±VL;X ZU 


97 Mnv 901 


600 


4.14 


600 


24.2 


RRT T97 

J_> ± VXj X Z 1 


9Q Mnv 901 
zy-iN u v-zuxu 


1600 


2.22 


1960 


27.6 


RRT T30 


9k Mnv 901 

ZO IN UV iUlU 


800 


6.39 


1160 


53.3 


RRT.T31 

J_X ± V 1 i X X 


1 8-Son-201 1 


2000 


6.18 


2360 


27.6 


RRT.T32 

!_} 1 \ 1 i X O ^1 


20-Sen-201 1 


2000 


5.76 


2360 


25 


RRT T33 


1 Q ^on 901 1 

X £7 UCU ZiUl X 


1600 


2.93 


1960 


20.9 


J_> ± \,Xj lOJ 


91 Son 901 1 

Z 1 UCU iUl X 


1600 


5.86 


1960 


21.6 


RRT T37 


1 8 Son 901 1 


2000 


4.84 


2360 


26 


RRT T38 


97 Nnv 901 


240 


4.14 


600 


31.9 


RRT TV) 

J_> ± VXj IOcj 


1 8 Son 901 1 


1600 


7.77 


1960 


29.8 


RRT T49 

X_> ± \,Xj X 1i 


90 ^on 901 1 


1600 


5.75 


1960 


23.6 


RRT T44 


30 Nnv 901 

OU IN U V-ZrUXU 


1600 


4.3 


1960 


28.6 


RRT T4 c i 


91 Spn 901 1 


1600 


3.76 


1960 


16.7 


RRT,T4fi 

i_} ± LXj X rr\J 


^O-Nnv-901 


1600 


3.84 


1960 


21.6 


RRT T48 


98 Nnv 901 

ZO IN U V ZiUXU 


1400 


2.93 


1760 


27.6 


RRT T^f) 
x_> ± \,xj x o u 


30 Nnv 901 

OU IN UV ZiUiU 


2400 


~1 


2760 


7.26 


RRT T^il 

J_> ± VXj X U X 


1 Q ^on 901 1 

X Ot[) ZrUX X 


2000 


4.59 


2360 


22.9 


RRT T c i9 

J_>±\,X^X uz 


97 Nnv 901 


1200 


6.55 


1200 


29.1 


RRT 

J_> ± i_Xj X u u 


9Q Nnv 901 


1600 


3.01 


1960 


22.8 


X_> ± \,Xj X O 1 


30 Nnv 901 

OU IN UV-ZiUlU 


2400 


3.93 


2760 


26.6 


J_> ± v,Xj IiJO 


97 Nnv 901 


1200 


5.09 


1200 


27.1 


RRT Tfifl 

j_>j.\x-jx uu 


9^ Ron 901 1 

1 l_- U ZiU X X 


1400 


4.11 


1760 


27.9 


RRT Tfi9 

J_> ± \,Xj X U Zi 


9Q Nnv 901 

Zii7 IN UV AUlU 


1600 


2.64 


1960 


26.8 


RRT TfI4 

J_>±\,XjX Urt 


97 Nnv 901 


1200 


5.26 


1200 


29.5 


RRRTfifi 

1_> ± LXj X \J\J 


27-Nnv-201 


800 


4.82 


800 


27.5 


RRT T^O 1 ^ 


Q 7 11T1 901 1 
O - J UXl-ZUX 1 


1 fiOO 
1UUU 


1 94 


1 QfiO 
xyuu 


1Q Q 


RRT T^Ofi 
n> jxxj x ouu 


7 Tun 901 1 


1600 


1.71 


1960 


25.7 


RRT T307 


1 8 Son 901 1 
X o otu ^Ul X 


2000 


7.84 


2360 


27.2 


RRRT31 1 

1_> 1 LXj 1 Ull 


1 Q-Sen-201 1 

X O kj ]X ^, U X X 


1600 


2.93 


1960 


16.2 


BRLT312 


20-Sep-2011 


2000 


4.38 


2360 


22.9 


BRLT313 


25-Feb-2011 


1600 


6.34 


1960 


42 


BRLT314 


7-Jun-2011 


1200 


1.71 


1560 


28.2 


BRLT315 


19-Sep-2011 


1600 


4.96 


1960 


25.3 


BRLT316 


21-Sep-2011 


2000 


4.41 


2360 


20.3 


BRLT317 


28-Nov-2010 


800 


13.7 


1160 


77.1 


BRLT318 


18-Sep-2011 


1600 


4.91 


1960 


21.7 


BRLT320 


20-Sep-2011 


1600 


4.49 


1960 


19.7 


BRLT321 


8-Jun-2011 


1600 


1.26 


1960 


16 


BRLT322 


18-Sep-2011 


1600 


5.48 


1960 


27.7 


BRLT323 


29-Nov-2010 


800 


3.83 


1160 


37.8 


BRLT325 


19-Sep-2011 


1600 


2.39 


1960 


18 


BRLT328 


29-Nov-2010 


1600 


4.23 


1960 


23.4 
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Figure 2. Left: Colour-colour diagram of candidate L and T dwarfs (crosses). Right: Colour-Magnitude diagram of candidate L and T 
dwarfs (crosses). Known L and T dwarfs (L0-T9) from Dwarfarchives.org are shown as filled red circles. Those with spectral types L3-T5 
dwarfs are shown as small blue diamonds. Photometry is on the UKIDSS and SDSS system. Objects with 2MASS photometry have been 
converted using the colour conversions of lWarren. Dye fc HambTvl <2006ft . 



Table l2l continued. Optical and near-infrared spectroscopic observations of sub-sample. 



BRLT330 


20-Sep-2011 


1600 


5.2 


1960 


20.6 


BRLT331 


19-Sep-2011 


2000 


5.45 


2360 


21.9 


BRLT332 


18-Sep-2011 


2000 


7.87 


2360 


27.7 


BRLT333 


28-Nov-2010 


800 


1.7 


1160 


21.5 


BRLT334 


27-Oct-2010 


240 


4.06 


600 


33.9 


BRLT335 


21-Scp-2011 


2000 


3.75 


2360 


24.5 


BRLT338 


27-Nov-2010 


800 


3.44 


800 


18.2 


BRLT343 


28-Nov-2010 


1200 


3.67 


1560 


28 


BRLT344 


18-Sep-2011 


2000 


4.17 


2360 


26.8 



*S/N is per spectral pixel, after binning. 



to the templates using the spectral typing schem e from 
iKirkpatrick et al.l jl999f) and iBurgasser et alj |2006j ), for L 
and T dwarfs, respectively. The three best fits were also vi- 
sually inspected to make sure that they were good fits to the 
real spectra. The spectral types, as well as their associated 
errors are shown in Table. [4] and the spectra can be seen in 
ascending spectral type in Fig. [3] - [9] 

5.2 Unresolved binarity 

The L-T transition region is thought to be hindered with 
a substantial number of unresolved binaries (up to 45%; 



iMaxted fc Jeffriesll2005h . such that in order to measure a re- 
liable formation history we must distinguish between these 
binaries. With the wavelength coverage of our X-shooter 
data we should be able to differentiate between binaries 
where the components are separated by more than 1 spectral 
type, leaving only the equal spectral-type binaries, which 
should not affect the overall shape of the T e g distribution, 
since the equal T G ff bina ry fraction should b e uniform over 
this spectral type range l|Geiffler et al.ll201ll ). 

We select potential binary cand idates using the spectra l 
indices and the criteria defined in IBurgasser eraD (|2010l ) 
who use index-index and index-spectral type combinations 
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Table 3. Spectral templates from the SpeX Prism Spectral Library used for spectral typing. 



Object Name 


2MASS Designation 


Optical SpT 


NIR SpT 


Resolution 


Reference 


VB 8 


T1 fi^^^Q 0893401 


M7 V 


M7 


120 


I 


V 1 J 1U 


T1 Q1 fi r »7fi9-i-0 c »0Q091 


M8 V 


M8 


120 


2 


LTTS 9994. 


T1 4984393-U331 03Q1 


M9 V 


M9 


120 


3 


9TvTA fc ^W T91 ^OzLAfi 084^90 


T91 3044K4 084^90^ 


LI 5 


LI 


120 




Kelu-1 


T1 30^401 Q 9 c i410^Q 


L2 


L2 


120 


5 


91VTA t ^ t lW T1 ^0^444-1 '391 Ofi 


T1 r ^nfi r >441 4-1 391 OflO 

lOUUJ'illy IOZIUUU 


L3 


L3 


120 


Q 


9A/TA f ^ c l 191^804^7 1 ^OOQS 


191^804^7 1 ^OOQK 


L4 


L i 


120 


4 


OIJOO J UOOdUU. lUTiydOU'i.'i 


T083 c i0fi99-Ul Q^O^O 




L5 


120 


7 


2MASSI J1010148-040649 


J10101480-0406499 


L6 


L6 


120 


8 


2MASSI J0103320+193536 


J01033203+1935361 


L6 


L7 


120 


9 


2MASSW J1632291+190441 


J16322911+1904407 


L8 


L8 


75 


5 


DENIS-P J0255-4700 


J02550357-4700509 


L8 


L9 


120 


10 


SDSS J120747.17+024424.8 


J12074717+0244249 


L8 


TO 


120 


11 


SDSS J015141. 69+124429.6 


J01514155+1244300 




IT 


120 


12 


SDSSp J125453.90-012247.4 


J12545393-0122474 


T2 


T2 


120 


12 


2MASS J12095613-1004008 


J12095613-1004008 


T3.5 


T3 


120 


12 


2MASSI J2254188+312349 


J22541892+3123498 




T4 


120 


12 


2MASS J15031961+2525196 


J15031961+2525196 


T6 


T5 


120 


12 


SDSSp J162414.37+002915.6 


J16241436+0029158 




T6 


120 


13 


2MASSI J0727182+171001 


J07271824+1710012 


T8 


T7 


120 


13 


2MASSI J0415195-093506 


J04151954-0935066 


T8 


T8 


120 


12 



1: Burgasser et al. (20081 . 2: Burgasser (20041 . 3: Burgasser fc McElwain (20061 . 4: Kirkpatrick et al. (20101 . 

5: Burgasser (20071 , 6: Burgasser et al. (20071 . 7: Chiu et al. (20061 . 8: Reid et al. (20061 . 

9: Cruz et al. (20041 , 10: Burgasser et al. (20101 . 11: Looper, Kirkpatrick fc Burgasser (20071 . 

12: Burgasser (20041 , 13: Burgasser et al. (20061 



Table 5. Binary indices as denned b y iBur gasser et al 

except for //-dip, which is defined by iBurgasser et al. r<i2oioh . 



Table 6. Inflect i on po ints used for analysis of unresolved binarity 
IBurgasser et all <201CH . 



Index 


Numerator 


Denominator 


Feature 




Range 


Range 




H2O-J 


1.14-1.165 


1.26-1.285 


1.15 lira H2O 


H2O-F 


1.48-1.52 


1.56-1.60 


1.4 jim H2O 


H 2 0-,R" 


1.975-1.995 


2.08-2.10 


1.9 Jim H 2 


CH4-J 


1.315-1.34 


1.26-1.285 


1.32 fj,m CH 4 


CH 4 -# 


1.635-1.675 


1.56-1.60 


1.65 ^m CH 4 




2.215-2.255 


2.08-2.12 


2.2 fim CH 4 


K/J 


2.060-2.10 


1.25-1.29 


J-K colour 


H-dip 


1.61-1.64 


1.56-1.59 + 1.66-1.69 


1.65 (an CH 4 



to define six different criteria to segregate possible unre- 
solved L-T transition binaries, using known double objects 
as reference. These spectral indices take into account the 
flux ratio of the prominent molecular absorption bands of 
water and methane over the flux emitted in the J, H and 
K bands. Objects that match two criteria are called weak 
candidates while objects that match at least three of them 
are called strong candidates. For clarity, all criteria used to 
classify these dwarfs are re-stated in Table. [5] and [6] We 
show our observed sample and the binary selection criteria 
applied in the different index-index and index-spectral type 
in Fig. 1101 In this way we identified six strong candidates and 
six weak candidates, which are highlighted as diamonds and 
asterisks, respectively. 

The 12 unresolved binary candidates were investigated 
further using synthetic binary template fitting. The syn- 
thetic binaries were created by taking the standard spectral 
type templates available in the SpeX-Prism library. These 



Abscissa 


Ordinate 


Inflection Points 

(x,y) 


H2O-J 


U2O-K 


(0.325,0.5), (0.65,0.7) 


CH 4 -# 


CH4O-.fi: 


(0.6,0.35), (1,0.775) 


CH 4 -# 


K/J 


(0.65,0.25), (1,0.375) 


H2O-F 


H-dip 


(0.5,0.49), (0.875,0.49) 


SpT 


H 2 0-J/H 2 0-i/ 


(L8.5,0.925),(T1. 5,0.925), (T3. 5,0.85) 


SpT 


H 2 0-J/CH 4 -if 


(L8. 5,0.625), (T4.5,0.825) 



were then normalized to 1 at 1.28jim and scaled the spectra 
to an absolute flux level, using t he absolute J magnit ude - 
spectral type relation derived in lMarocco et afl (|2010l ). Fi- 
nally all the templates were combined to make our own set 
of synthetic unresolved binaries. We then performed a chi- 
squared fit of our spectra to the templates. The results of 
the fitting are presented in Table. 15.21 with the best-fitting 
combined template over-plotted on the spectra of the 12 bi- 
nary candidates shown in Fig. 1111 and 1121 For each object 
we list the best fit single object template and the associated 
chi-squared, the best fit combined template and its associ- 
ated chi-squared fit value. The results of the two fits were 
compared using a one-sided F test to assess their statistical 
significance. If the ratio of the two chi-squared fits (rj) is 
greater than the critical value (1.15), this represents a 99% 
significance that the combined template fit is better than 
the standard template alone. The last column of Table 15.21 
presents the results of the F test, where it can be seen that 
six of the 12 candidates give statistically better fits using 
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Table 4. Spectral details of observed targets. 



ID Spectral type H 2 Q-J H 2 0-H H a O-K CH 4 -J CH 4 -H CH 4 -K K/J H-dip Notes 



BRLT1 


L9.0 


± 


0.5 





72 


0. 


.68 





,86 


0. 


75 


1 


.00 


0.76 





.60 





.49 




BRLT2 


L1.0 


± 


1.0 





.94 


0, 


76 


1 


.13 


0. 


80 


1 


.09 


1.06 





.40 





,46 




BRLT3 


L9.0 


± 


1.0 





.75 





,69 





.87 





.73 


1 


.04 


0.90 





.58 





.50 




BRLT6 


L3.0 


± 


1.0 


0. 


.85 





,82 


1 


.04 


0. 


.83 


1 


.13 


1.11 





,45 





.49 




BRLT7 


M8.0 


± 


1.0 





.95 





.87 


1 


.15 





.86 





.98 


0.91 





.36 





.50 




BRLT8 


L8.5 


± 


0.5 





71 


0, 


75 





.87 


0, 


,77 


1 


.03 


0.89 





.67 





,49 




BRLT9 


L1.0 


± 


1.0 





.95 





84 





,98 





,86 


1 


.03 


0.95 





.36 





,48 




BRLT10 


L9.0 


± 


0.5 





.67 


0, 


72 





,79 


0, 


76 


1 


.05 


0.85 





.64 





,52 




BRLT12 


L3.0 


± 


1.0 





,75 


0, 


76 





,95 


0, 


,74 


1 


.05 


0.99 





.45 





,47 




BRLT14 


LO.O 


± 


0.5 





97 


0, 


84 


1 


,13 


0, 


86 


1 


.05 


1.05 





.37 





,48 




BRLT15 


T2.0 


± 


2.0 





49 





.63 





,81 





.57 





.92 


0.75 





,52 





.47 


SC 


BRLT16 


L3.5 


± 


0.5 





.90 





,75 





,96 





.83 





.93 


0.95 





.42 





.45 


WC 


BRLT18 


LO.O 


± 


1.0 


1 


03 


0, 


91 


1 


,02 


0. 


94 


1 


.03 


0.95 





.39 





,49 




BRLT20 


L1.0 


± 


1.0 





.87 


0, 


76 


1 


,00 


0. 


,76 





.99 


0.91 





.32 





,47 




BRLT21 


L3.5 


± 


0.5 


0. 


83 


0, 


,75 





.95 





,80 


1 


.00 


0.98 





,45 





,47 




BRLT22 


M9.5 


± 


0.5 


1 


.01 


0, 


,90 


1 


.16 





.87 


1 


.12 


0.86 





.32 





.51 




BRLT24 


L3.5 


± 


0.5 





.80 


0, 


74 





,96 





,80 


1 


.00 


0.97 





.47 





.46 




BRLT26 


L5.5 


± 


0.5 


0. 


80 





,72 





,91 


0. 


,81 


1 


.07 


0.96 





.53 





,48 




BRLT27 


T1.0 


± 


0.5 





62 





,64 





,78 





.63 


1 


.01 


0.59 





.40 





.50 




BRLT30 


L5.0 


± 


0.5 





.74 


0, 


,75 





,95 


0, 


.79 


1 


.08 


1.01 





.51 





,48 




BRLT31 


L4.0 


± 


1.0 





,74 


0, 


69 





,92 


0, 


,77 


1 


.01 


0.85 





.51 





.50 




BRLT32 


LI. 5 


± 


0.5 





.85 


0, 


79 


1 


,08 


0, 


83 


1 


,04 


1.05 





.40 





.48 




BRLT33 


L3.5 


± 


0.5 





,82 





,73 





,90 





,80 





.94 


0.91 





,45 





.47 


wc 


BRLT35 


M9.5 


± 


0.5 


1 


.01 


0, 


,85 


1 


.00 





.93 


1 


.04 


1.00 





.41 





,48 




BRLT37 


L5.0 


± 


0.5 





71 


0, 


71 





.82 


0, 


,76 


1 


.09 


0.93 





.53 





,50 




BRLT38 


T1.0 


± 


0.5 





65 


0, 


65 





,77 


0, 


68 





.95 


0.67 





.40 





.51 




BRLT39 


L5.0 


± 


1.0 





,79 


0. 


.75 





,93 





,77 


1 


.04 


0.95 





.46 





.50 




BRLT42 


M9.0 


± 


0.5 


1 


,04 


0. 


,89 


1 


.10 





.95 





.94 


0.94 





.45 





,47 


WC 


BRLT44 


L5.0 


± 


1.0 





72 


0, 


,75 





.83 


0, 


,74 


1 


.00 


0.81 





.52 





.49 




BRLT45 


T1.0 


± 


0.5 





57 


0, 


61 





.71 


0, 


64 





.88 


0.57 





.42 





.47 


SC 


BRLT46 


L0.5 


± 


0.5 





,91 





,76 


1 


.13 


0, 


,80 


1 


.12 


1.07 





.30 





,48 




BRLT48 


L4.5 


± 


0.5 





,77 





,77 


1 


.04 





,81 


1 


.19 


1.08 





.55 





,49 




BRLT50 


T7.0 


± 


0.5 





.16 





,33 





,84 


0, 


,28 





.25 


-0.0 





.06 





.20 




BRLT51 


L3.0 


± 


1.0 





.80 


0, 


,78 





,86 


0, 


83 


1 


.06 


0.92 





.51 





,47 




BRLT52 


L5.5 


± 


0.5 





75 


0, 


68 





.94 


0, 


,77 


1 


.05 


0.94 





.57 





.50 




BRLT56 


LI. 5 


± 


1.0 





,91 





86 





.93 


0. 


,85 


1 


.01 


0.96 





.14 





.49 




BRLT57 


LO.O 


± 


1.0 





,97 


0, 


,86 


1. 


.08 





.84 


1 


.09 


0.99 





.37 





.50 




BRLT58 


L5.5 


± 


1.0 





76 


0, 


72 





,84 


0, 


,74 


1 


.05 


1.01 





.41 





.49 




BRLT60 


L1.0 


± 


1.0 





,83 


0, 


79 





,99 


0, 


,84 





.99 


0.94 





.41 





.46 


wc 


BRLT62 


L5.0 


± 


1.0 





,83 





,76 





.97 





,86 


1 


.11 


1.04 





.63 





.49 




BRLT64 


L4.0 


± 


0.5 





,76 





,73 





.92 





,83 


1 


.05 


1.01 





.50 





,48 




BRLT66 


L5.0 


± 


0.5 





74 


0, 


69 





,92 


0, 


79 


1 


.08 


0.95 





.52 





,48 




BRLT305 


L5.5 


± 


1.0 





.93 


0, 


79 





.95 


0, 


91 





.97 


0.93 





.61 





.47 


wc 


BRLT306 


L4.5 


± 


1.0 


0. 


81 


0, 


,80 





.95 





,75 


1 


.01 


0.92 





,39 





.48 




BRLT307 


L1.0 


± 


0.5 





,92 


0, 


,89 


1 


17 


0, 


,93 


1 


.07 


0.93 





.47 





,48 




BRLT311 


T3.0 


± 


0.5 





42 


0, 


56 





.78 


0, 


.55 





.70 


0.57 





.29 





,44 


SC 


BRLT312 


TO.O 


± 


0.5 


0. 


62 


0, 


67 





.87 


0, 


66 





.92 


0.90 





.44 





.42 


SC 


BRLT313 


L3.5 


± 


0.5 





,82 


0, 


,74 





.94 





,88 


1 


.04 


0.99 





.53 





.47 




BRLT314 


L7.5 


± 


0.5 





,72 





,71 





,88 





,80 


1 


.05 


0.95 





.69 





,49 




BRLT315 


L1.0 


± 


1.0 





.87 


0, 


82 





,92 


0, 


,78 


1 


.04 


1.01 





.35 





,47 




BRLT316 


L1.0 


± 


0.5 





.87 


0, 


87 





,91 


0, 


,81 


1 


.13 


0.98 





.35 





.49 




BRLT317 


L1.0 


± 


1.0 


0. 


,79 





,71 





.92 





,70 


1 


.09 


0.95 





.31 





,49 




BRLT318 


L1.0 


± 


0.5 





83 


0, 


,88 





.97 





.85 


1 


.13 


0.95 





.42 





.49 




BRLT320 


L1.0 


± 


0.5 





.89 


0, 


79 





.90 


0, 


,85 


1 


.10 


0.89 





.42 





.50 




BRLT321 


T4.0 


± 


0.5 





.37 


0, 


42 





.57 


0, 


43 





.60 


0.25 





.24 





.37 




BRLT322 


L5.0 


± 


0.5 





,71 


0, 


,70 





.90 





,72 


1 


.07 


0.90 





.55 





.52 




BRLT323 


TO.O 


± 


1.0 





.76 





,71 





.87 


0, 


,76 


1 


.04 


1.04 





.42 





.49 




BRLT325 


T2.0 


± 


1.0 





,53 


0, 


.47 





.63 


0, 


55 





.89 


0.49 





.27 





,48 




BRLT328 


L0.5 


± 


1.0 





.89 


0, 


78 





,97 


0. 


,82 


1 


.09 


1.05 





.36 





,50 
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Table l4l continued. Spectral details of observed targets. 



ID Spectral type H 2 Q-J H 2 Q-H H 2 Q-K CH 4 -J CH 4 -H CH 4 -K K/J H-dip Notes 



BRLT330 


L1.0 ± 1.0 


0.81 


0.77 


0.89 


0.73 


1.20 


0.95 


0.35 


0.48 




BRLT331 


L1.0 ± 2.0 


0.84 


0.71 


1.06 


0.74 


1.05 


0.93 


0.30 


0.49 




BRLT332 


L2.0 ± 2.0 


0.87 


0.85 


1.00 


0.75 


1.10 


0.96 


0.34 


0.49 




BRLT333 


T2.0 ± 0.5 


0.52 


0.57 


0.77 


0.59 


0.98 


0.72 


0.38 


0.52 


SC 


BRLT334 


L3.5 ± 0.5 


0.83 


0.74 


0.91 


0.82 


1.04 


0.95 


0.48 


0.48 




BRLT335 


L4.0 ± 1.0 


0.81 


0.74 


0.95 


0.79 


0.99 


0.92 


0.42 


0.47 


wc 


BRLT338 


L1.0 ± 1.0 


0.88 


0.73 


1.01 


0.81 


1.09 


1.12 


0.34 


0.49 




BRLT343 


L9.0 ± 1.0 


0.69 


0.70 


0.88 


0.75 


1.05 


0.92 


0.60 


0.49 




BRLT344 


TO.O ± 1.0 


0.54 


0.68 


0.96 


0.68 


0.94 


0.89 


0.54 


0.44 


SC 



SC=Strong binary candidate, WC= Weak binary candidate 



♦ 

♦ *l 



♦ ♦ 



♦ . 



Strong candidates ♦ 
Weak candidates * 



0.4 0.6 0.8 

HO-J 
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Strong candidates ♦ 
Weak candidates * 



0.5 0.6 0.7 0.8 0.9 1.0 1.1 
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0.5 
0.4 
0.3 
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0.42 
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i Weak candidates * - 


0.40 
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0.6 
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NIR Spectral Type 



•j. 



r 
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NIR Spectral Type 



Figure 10. Index-Index plots showing the selection criteria applied for unresolved binary candidate segregation. Strong and weak 
candidates are marked as diamonds and crosses respectively. 



the combined templates, such that they are the strongest 
binary candidates. It should b e noted that the resu lts of 
this fitting are not conclusive l|Burgasser et alj |2010), and 
the real nature of these objects must be investigated fur- 
ther using adaptive optics imaging or spectro-astrometry in 
order to confirm or not their true binary nature. 



6 CONSTRAINING THE GALACTIC BROWN 
DWARF FORMATION HISTORY 

In order to compare the findings from our first sub-sample of 
76 mid L to mid T dwarfs that occupy a complete area of 495 
sq degrees of sky (RA=22 to 4 hr and DEC=-2 to 16 deg) 
down to a magnitude limit of J=18.1, with the results from 



simulations of differing birth rates, we constructed a space 
density vs spectral type/T G ff histogram. We firstly calculated 
the space density of our targets by converting their spectral 
types into effective temperatu res using the T e g— NI R spec- 
tral type relation presented in lStephens et aL I (120091 , equa- 
tion 3). Using a bin sizes of 1700-1450K, 1150-1300K and 
1300-1450K (corresponding to bins in spectral type of L4-L6, 
L7-T0 and T1-T4), we calculated for each bin the maximum 
distance at which an object could have bee n selected using 
the M j— NIR spectral type relation from iMarocco et al.l 
(2010). With this distance limit we then calculated the vol- 
ume sampled by each T e a bin. The derived space densities 
were then corrected for Malmqui st and Edding t on bi as fol- 
lowing the approach described in lPinfieid et al.l(|2008l ). The 
Eddington bias is caused by the photometric uncertainties 
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BRLT31 1 

2MASSJ1209b613-100400B (T3) 
T1 .0+T5.0 

SDSSJ015141. 69+124429.6 (T1.0) 
2MASSJ15031961+2525'96 (T5.0) 



3 r 




1 fwyw" 



Wavelength (jU-m) 



BRLT333 

SDSSpJ125453. 90-012247, 4 (T2) 
T0.0+T2.0 

SDSSJ120747. 17 + 024424. 8 (TO.O) 
SDSSpJ125453. 90-012247. 4 (T2.0) 




Wavelength (/Am) 



BRLT312 

S0SSJ120747. 17 + 024424. 8 (TO) 
L5.0 + T4.0 

SDSSJ083506. 16+195304. 4 (L5.0) 
2MASSIJ2254188 + 312349 (T4.0) 




Wavelength (/xm) 



BRLT344 

SDSSJ120747. 17 + 024424, 8 (TO) 
L7.0 + T7.0 

2MASSIJ0103320+'93536 (L7.0) 
2MASSIJ0727182 + 171Q01 (T7.0) 




Wavelength (/2m) 



Figure 11. "Strong" unresolved binary candidates over-plotted with combined spectral templates. 



on the magnitudes of objects near our cut (i.e. J < 18.1). 
However, since the magnitude cut we impose is bright (it 
corresponds to a ~ 12a detection in the UKIDSS LAS), 
the uncertainties at the J=18.1 limit are typically less than 
a=0.05 and therefore the Eddington bias correction is less 
than 1%. This is negligible compared to the other sources 
of uncertainty. We estimated the Malmquist bias correction 
considering the mean scatter of the sample of known L and 
T dwarfs around the adopted Mj— NIR spectral type rela- 
tion. This represents an increase in the volume sampled of 
22%. 

6.1 Completeness 

In order to calculate the completeness of our sample we com- 
pared it to a control sample of known L and T dwarfs taken 
from DwarfArchives.org, for a magnitude limit of J ^ 16, 



removing any objects that are known to be members of 
unresolved binary systems. The control sample was cross 
matched with the UKIDSS LAS and SDSS in order to ob- 
tain photometry on the same colour system as our selection 
criteria. We use DR9 of the UKIDSS LAS and DR8 of the 
SDSS in order to get a sufficiently large control sample. We 
imposed the same set of colour cuts, as described in to 
reveal the level of completeness of our sample selection. We 
retain all of the L4 dwarfs from the control sample, but only 
some of the L0-L3 dwarfs when imposing our colour cuts, 
as such we find that our sample selection is complete for L4 
spectral types and later. 

We also consider how many objects we may lose from 
our selection due to photometric scattering of colours. We 
calculate that for L0-L3 types we would expect to lose 3.7 
dwarfs, this corresponds to a completeness level of 85% for 
L0-L3 dwarfs due to photometric scatter, although we note 
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BRLT16 

2MASSWJ1506544+132106 (L3) 
L3.0+T3.0 

2MASSWJ1506544+132106 (L3.0) 
2MASSJ12095613-1004008 (T3.0) 



T3 3 
CD 

E 2 
o 




Waveleng;h (/xrn) 



£ 2 W* 8 ^ 



BRLT33 

2lv'ASSWJ 1506544+ 132 106 (L3) 
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Figure 12. " Weak" unresolved binary candidates over-plotted with combined spectral templates. 



that we are significantly incomplete from our colour selection 
for this spectral range. For L4-L6 dwarfs we would expect 
to lose 2.33 dwarfs, being 88% complete. The L7-T0 range 
would lose 0.55 dwarfs, corresponding to a 94% complete- 
ness; and for T1-T4 we would expect to lose 0.05 dwarfs, 
corresponding to a completeness of 99%. 

We combine the completeness for colour selection, 
colour scatter as well as the Malmquist and Eddington biases 
to infer a 90% completeness for the T e g =1450-1700K and 
97% for the T eff = 1300-1450K and 97% for the T cft = 1150- 
1300K range in our sample. 

6.2 Correction for unresolved binarity 

We also corrected our results for the presence of binaries by 
firstly removing objects identified as possible binaries ( £|5. 2|) 
for which our spectral deconvolution gives a statistically bet- 



ter fit. We then applied a further correction to take into 
account the presence of equal spectral type binaries, which 
would fall beyond our J < 18.1 limit if they were single 
objects, using the definit i on of "observed binary fraction" 
given bv lBurgasser et all (^OOSt ). such that: 



N B 



7 + 1 
BF — 1 



(1) 



where Nb and N m are the observed binaries and the 
total number of objects respectively, BF is the "true" binary 
fraction, and 7 is the fractional increase in volume due to 
inclusion of binaries in the sample. The number of binaries 
that fall within our magnitude limit (Nd) is: 



N B 



7' 



(2) 
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Figure 3. X-shooter spectra of L and T dwarfs from birth rate 
sample, in ascending spectral type. 



Figure 4. X-shooter spectra of L and T dwarfs from birth rate 
sample, in ascending spectral type. 



Table 7. Results of spectral fitting of binary candidates using 
combined templates. 



Target 


Single template 


Combined template 


F-test 


ID 


best fit (x 2 ) 


best fit (x 2 ) 


V 




Strong 


candidates 




BRLT15 


T2.0 (14.8341) 


L8.0+T7.0 (3.6926) 


4.02 


BRLT45 


T1.0 (6.2583) 


T1.0+T3.0 (6.7607) 


0.92 


BRLT311 


T3.0 (9.3321) 


T1.0+T5.0 (7.2732) 


1.28 


BRLT312 


TO.O (6.5036) 


L5.0+T4.0 (5.3935) 


1.20 


BRLT333 


T2.0 (3.5947) 


T0.0+T2.0 (3.7963) 


0.95 


BRLT344 


TO.O (5.8038) 


L7.0+T7.0 (3.5877) 


1.62 


Weak candidates 


BRLT16 


L3.5 (3.7444) 


L3.0+T3.0 (3.6042) 


1.04 


BRLT33 


L3.5 (6.0794) 


L3.0+T4.0 (4.4195) 


1.37 


BRLT42 


M9.0 (5.2356) 


M9.0+L1.0 (3.3914) 


1.54 


BRLT60 


L1.0 (2.2994) 


L1.0+T3.0 (2.1760) 


1.06 


BRLT305 


L5.5 (7.0290) 


L5.0+L6.0 (8.9120) 


0.79 


BRLT335 


L4.0 (4.6314) 


L3.0+T3.0 (4.2767) 


1.08 



Therefore the fraction of objects to be excluded from 
our sample (f exc i) is 

1 



7" 



_ Nb_No = 
Uxcl ~ N m N B ~ 7 + 1/BF - 1 



(3) 



For equal spectral type binaries 7 = 2v2. For the "true" 
binary fraction we assumed two extremes valu e s take n 
from the literature. Following iBurningham et ail (|2010al ). 
for the lower limit w e considered the values obtained by 
iBurgasser et alj (|2003l ) who estimated a BF of 5-24% using 
A O imaging; for the upper limit we used the values given 
bv lMaxted fc Jeffries! l|2005l ), who estimated a BF of 32-45% 
via radial velocity monitoring. 

6.3 Space densities 

We combine our calculated space densities, taking into ac- 
count our completeness and contribution from un resolved 
binaries, with those fromlBurningh am et al.l (l2010aT) for late- 
type T dwarfs. We choose to use the Burningham et all 
(2010a) space den sities in favour of others available in 
the literature (e. g. ICruz et all l200Tt iMetchev et all 120081 . 
iRevle et al.l I2010L iKirkpatrick et alj 12012^ " as this is the 
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Figure 5. X-shooter spectra of L and T dwarfs from birth rate 
sample, in ascending spectral type. 

only space density that can be used in direct compari- 
son to our own sample, as we use the same Mj-spectral 
type relations as well as apply the same binary fractions. 
It is also probes d own to a magnitu de limit comparable 
to our sample. The I Cruz et all (2007) space densities only 
probe out to the 2MASS limit (~ J = 16) and covers 
only M9-L8 dwarfs, and not the full range of spectral types 
we are probing. The iMetchev et al.l (|200ct ) sample provide 
space densities obtained from 2MASS, combined with SDSS, 
and as such does not probe as deeply as our sample, or 
completely into the late T dwa rf spectral type as that of 
iBurningham et all l|2010ah . The lRevle et al.l l|201Cl ) derived 
space densities could also be comparable but for their very 
late T type densities (T6-T8), they only have a handful of 
obj ects and do not make any c orrection for binarity. As such 
the IBurningham et al] l|2010ah space densities act as a com- 
plementary extension to our mid L-mid T derived space den- 
sities. Our calculated space densities are shown in Table[8l 
with comparison to the above mentioned published space 
densities for L and T dwarfs. 

Our space densities are also shown as a function of both 
spectral type in Fig. 1131 and as a function of T e a in Fig. 1141 
where we show two points for each of our T e g /spectral type 



E 
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Wavelength (fim) 



Figure 6. X-shooter spectra of L and T dwarfs from birth rate 
sample, in ascending spectral type. 

bin, the upper most points represent the density correspond- 
ing to a binary fraction of just 5%, while the lower points 
are those for a binary fraction of 45%. Also overplotted are 
space densities published for L and T dwarfs. In addition 
we show (overplotted) the results of numerical simulations 
computed assuming different IMF and birth rates. De t ails o f 
the simulations are presented in iDeacon fc Hamblvl (|2006t ) 
and are briefly summarized here. 

We assume an exponential IMF in the form: 

ip(M) oc M~ a (pc^Mo -1 ). (4) 

where ^ is the number of objects per unit volume in a 
given mass interval. We also assumed an exponential birth 
rate of the form: 

b(t) oc e" /3t (5) 

where t is in Gyr and /3 is the inverse of the scale time 
t (in Gyr, since the galaxy was formed). Each simulated 
object was assigned an age based on the birth rate and a 
mass based on the IMF, giving a final creation function C 
given by the equation: 

C(M,t) = *(M)^ (6) 
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Figure 7. X-shooter spectra of L and T dwarfs from birth rate 
sample, in ascending spectral type. 



Figure 8. X-shooter spectra of L and T dwarfs from birth rate 
sample, in ascending spectral type. 



where To is the age of the Galaxy. C is therefore the 
number of objects created per unit time per unit mass. 

The evolution of each object and its parameters (i.e. T e s 
and absolute magni tudes) were calculat ed using the evolu- 
tionary models from lBaraffe et all l|200&t h We note that any 
any model dependent systematics would be introduced, but 
that these should not effect the overall trend. The num- 
ber densities obtained for each temperature bin were finally 
normalized t o 0.0024 pc~ 3 in the 0.1—0.09 Ma m ass range, 
according to lDeacon. Nelemans fc Hambfvl l|200Sf ). 

We consider five different values of f3 (-0.2,- 
0.1,0.0, +0.1, +0.2 corresponding to r=-5,-10,oo,+10,+5 Gyr 
respectively) and five values of a (0.0,-0. 5,-1.0-, 1.5,-2.0). The 
results obtained for a = 0.0,-1.0,-2.0 and /3 = 0.0,0.2,0.5 are 
shown in Figure 1141 where different colours represent dif- 
ferent values of a and different line styles represent different 
values of /3. 

A first look at our measured mass function shows 
that the more extreme birth rates (e.g. a halo type 
form) can likely be excluded as the true form of the 
birth rate. Our space densities are in general agreement 
with the previous Ba y esian analysis of 2MASS L and T 
dwarfs by I Allen et ail l|l995h . and do not differ drastically 



(within our uncertainties) from those mas s functions previ- 
ously measured and discussed earlier (i .e. Cruz et al. 2007; 
jMetchev etaH 120081: iRevle et ail l2010l : iBurningham et all 
l2010al ; lKir"kpatrick et al.ll2012l ). The differences between our 
derived densities and those previously published can be ac- 
counted for in large due to the use of different Mj-SpT con- 
versions by the different groups. 

Our space densities do however suggestive of a pref- 
erence for an a < —1.0, which is consistent w it h the 
stud ies of late T dwar f s alone from iPinfield et alj (|2008l ) 
and IBurningham et all l)2010al ) , where they conservatively 
find a < 0. Indeed all three samples are indicative of an 
— 1 < a < 0. It can be seen however that it is not cur- 
rently possible to place robust constraints on the birth rate 
with this sub-sample. One of the largest sources of uncer- 
tainties is the binary fraction. This could be resolved with 
the follow-up of our unresolved binary candidates, by either 
adaptive optics imaging or radial velocity. In addition im- 
proved numbers from the completion of follow-up of our full 
sample would allow us to place better constraints on the 
birth rate and would allow us to rule out at least j3 = or 
/3 — 0.5 as the uncertainties are reduced. With spectroscopic 
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Table 8. Calculated space densities. 



t~i r>TOTT>Tl ft 1 


Sp. 


Type 




1\L CL111HJ Ltlo L 1 n ] 1 1 1 1 1 M 1 .1 1 1 1 


Binary 


Space density 




range 


(no correction) 


bias correction (%) 


correction (%) (xlO 3 pc 3 ) 








(xl0~ 3 pc~ 3 ) 






This paper 


L4 - 


L6.5 


0.345 ± 0.061 


22 


3-45 


0.176 ± 0.031 - 0.295 ± 0.052 




L7 - 


T0.5 


0.293 ± 0.050 


22 


3-45 


0.140 ± 0.024 - 0.235 ± 0.040 




Tl - 


T4.5 


0.235 ± 0.088 


22 


3-45 


0.106 ± 0.040 — 0.178 ± 0.067 


Burningham et al. (2010a) 


T6 - 


T6.5 


_ 


12-16 


3-45 


0.30 ±0.2 -0.59 ±0.39 




T7 - 


T7.5 


- 


12-16 


3-45 


0.40 ± 0.28 - 0.79 ± 0.55 




T8 - 


T8.5 


- 


12-16 


3-45 


58 + 051 — 1 1 +1 




T9 








Qi-LOQ i _|_ e; 7 


Reyle et al. (2010) 


L5 - 


L9.5 


- 


- 


_ 


2.0+™ 




TO - 


T5.5 


- 


- 




1.4+° 3 

-0.2 




T6 


- T8 


_ 


_ 




5.3_ 2 ' 2 


Cruz et al. (2007) 


M7 - 


M9.5 








4.9 ±0.6 




LO 


- L3 








1.7±0.4 




L3.5 


- L8 








2.2 ±0.4 


Metchev et al. (2008) 


TO - 


T2.5 






~50 






T3 - 


T5.5 






~22 






T6 


- T8 






~14 




Kirkpatrick et al. (2012) 


T6 - 


T6.5 






30 


1.1 




T7 - 


T7.5 






30 


0.93 




T8 - 


T8.5 






30 


1.4 




T9 - 


T9.5 






30 


1.6 




YO - 


Y0.5 






30 


1.9 



follow-up of a full sample the errors could be reduced by up 
to 50%. 



7 SUMMARY 

We have identified 63 new brown dwarfs that lie in the L-T 
transition region, including the identification of 12 possible 
unresolved binary systems, fn order to confirm these as real 
binary systems additional adaptive optics/ radial velocity 
measurements are planned. We converted our 76 strong new 
sample onto a histogram of spatial density vs T e ff /SpT and 
compared them with simu lations of differing birth rates from 
iDeacon fc HambTvl l|2006l ) . This showed that the mid L-mid 
T range are in general agreement (within the e rrors) with the 
indications from the late T dwarfs alone from IPinfield et al] 
(|2008h and iBurningham et al] l|2010aT ). such that a < 0. 
Indeed both our L-T transition sample and the late T dwarf 
samples are suggestive of — 1 < a < 0. A better constraint on 
the binary fraction, as well as a larger sample size however, 
is required before we can place robust constraints on the 
form of the brown dwarf birth rate, other than to suggest 
that a halo form of the birth rate is extremely unlikely. 
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